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Abstract 

We update the distance priors by adopting Planck TT,TE,EE + lowP data released in 
2015, and our results impose at least 30% tighter constraints than those from Planck TT+ 
lowP. Combining the distance priors with the combination of supernova Union 2.1 compi¬ 
lation of 580 SNe (Union 2.1) and low redshift Baryon Acoustic Oscillation (BAO) data, 
we constrain the cosmological parameters in the freely binned dark energy (FBDE) and 
FBDE-|-f2fc models respectively, and hnd that the equations of state of dark energy in 
both models are consistent with w = —1. Furthermore, we show that the tension with 
the BAO data at z = 2.34 from Lya forest (LyaF) auto-correlation and Combined LyaF 
cannot be relaxed in the FBDE and FBDE-I-Dfc models. 
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1 Introduction 


The current cosmic acceleration DEI of the universe implies that there is a mysterious 
component inventory in the universe, namely the dark energy (DE). Superbly accurate 
data are necessary for understanding the nature of DE. Although the usual luminosity- 
distance method based on supernovae of type la (SNIa) can impose certain constraints on 
DE, there are limitations for SNIa when we consider the nonzero spatial curvature 
or an evolving dark energy model. Because the measurements of SNIa span a relatively 
narrow range of redshifts, some other cosmological observations are called for. To impose 
constraints on DE models that are beyond SNIa’s domain, the acoustic peak method is 
usually taken as a geometric complement. The acoustic peaks can be imprinted onto not 
only the late-time power spectrum of cosmic microwave background (CMB) anisotropy, 
but also that of the non-relativistic matter from which the baryon acoustic oscillation 
(BAO) can be measured. Although the acoustic features are thought as one of the most 
powerful ways to probe DE, there are also some limitations for them. One limitation is 
that the effects on the matter correlations are weak. The other is that the evolution of 
linear density perturbations, including perturbations of both traditional components and 
DE, should be took into account when we £t a certain DE model to CMB data through 
the full Boltzmann analysis isiiniii]- However, even though this prerequisite is ready, this 
method will be prohibitive because of the long computing time. 

Compounding the situation is that for several more imaginative DE models, we can’t 
set up the equations for linear density perturbations or we can set up these equations 
only in principle. DGP model [S], for example, is given hrst by an action. If one wants to 
get the CMB temperature power spectrum from scratch using it, one will get stuck with 
complications of setting up these equations [Siiiniiiiiiig and treating them numerically. 
Another hard case is that some DE models are based on phenomenological considerations, 
and then there is no enough information to calculate the density perturbations. Caedas- 
sian model [131 E]) example, with a modihed Friedmann equation gets its energy 
density perturbations through the fluid flow approach of Hawking [151 CB] • However, the 
results are very different when we choose a different gauge, which creates a problem of 
understanding perturbations in the current Universe. 

Fortunately, some pioneers have developed a substitution for full Boltzmann analysis 
of CMB anisotropy, which can not only incorporate as much empirical information as 
possible but also avoid the calculation of the evolution of of linear density perturbations. 
The fundamental principle is that using certain characteristic distance scales to summarize 
the CMB data, namely the shift parameter R [13 HE] that determines the amplitude of 
acoustic peaks in the power spectrum of CMB temperature anisotropy and the acoustic 
scale IA that determines the acoustic peak structure [IH]. These two distance scales from 
CMB are so-called distance priors. From WMAP [201 El] to Planck satellite [22] , almost 
every data release was accompanied by a corresponding set of distance priors. 


2 


In this paper, we will update the distance priors with the latest CMB data [23], in 
particular including high £ polarizations data, released by Planck collaborations recently, 
and make comparison to the other two sets [22l [2l| derived from Planck 2013 data and 
Planck 2015 data excluding high £ polarizations data. The rest of this paper is arranged 
as follows. We describe our methodology and present our results in Sec. 2. We use the 
distance priors derived in this paper to constrain the equation of state of dark energy in 
Sec. 3. The conclusion and discussion are given in Sec. 4. 


2 Reconstruct the distance priors from Planck 2015 
data 


2.1 Methodology 


Distance priors can be directly derived from public data, such as Planck or WMAP data 
which are gained through a full Boltzmann analysis of CMB data. In order to obtain 
the distance priors, one should assume a cosmological model hrst |2S|. It is necessary 
to show that the distance priors are effective observables. We will see that it is indeed 
the case later. Moreover, since we usually use distance priors to deal with the late-time 
expansion history of the universe, we also add the baryon density today klhh^ to the 
distance priors, which is useful to probe the late-time universe but not sensitive to the 
cosmological models. 

We derive the distance priors I a and R by following [20l [26] . The comoving scale of the 
hrst acoustic peak is well determined as Ap = rs^z^)/^, where is the comoving sound 
horizon which is given by 


rs{,z) 

3 


c 


U/ll+z) 


da 


0 Jo 


a?E{a)^^ 'd){l + 




31500(TcMB/2.7iP)-", Tcmb = 2.7255iP . 


( 1 ) 


Here z* is the redshift to the photo-decoupling surface, which is given by the CAMB pack¬ 
age. If we observe an angle 9a subtending the transverse comoving scale Ap, the angular 
diameter distance is given as Da{z^) = • It’s obvious that different cosmological 

models will give different Ap and Da{z^:), but (1 -1- z^) and 9a are observables. Therefore, 
we can easily construct an important derived parameter, i.e. the acoustic scale Ia given 
by 


Ia = 


1 


(1 -|- z^) 


7iDa{z^) 

Tsiz^) 


( 2 ) 
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which is an effective observable. Here we used the angular diameter distance 




'/^sinn 



( 3 ) 


where sinn(a;) = sin(a;), x, sinh(a;) for flk < 0, Qk = 0, flk > 0. Here E{z) is given by 
E{z) = H{z)/Ho, i.e. 


E(z)^ 


r2j.(l + z)^ + + fl, 


de 


Pde{0) 


where is the present fractional radiation density 




1 ff -^eq 


Zeq = 2.5 X (Tcms/ 2 . 7 K) 


-4 


(4) 

(5) 


For ACDM and wCDM models, Pde{z)/Pdei^) equals 1 and (1 + respectively. 

Similarly, we can work out the other important derived parameter, i.e. ‘shift parameter’ 
R, which is also an effective observable. At low redshift, a comoving scale extending 
along the line of sight spanning dz is given by A = cdz/H{z). It’s obvious that different 
cosmological models will give different A and H{z), but dz is an observable. Therefore, 
we get a variable dz = AffM jf apply it to the decoupling epoch, we can construct 


R{z.) = (l + ^.)gAWg(^.) , (6) 

c 

Actually, what we usually use is its another famous version, namely the ‘shift parameter’ 

R dZlIISl 


R{Z,) = + 

c 

That is to say, R is just the traditional construction, but not the only one. 

In summary, Ia characterizes the CMB temperature power spectrum in transverse 
direction, and different I a gives different distribution of peaks and troughs in the spectrum; 
R characterizes the CMB temperature power spectrum in line-of-sight direction, and 
different R will magnify or reduce the amplitude of the acoustic peaks. 


2.2 Results 

Although the distance priors have been derived from Planck TT + lowP data in [22] , the 
constraints on the distance priors are weaker than those given by Mukherjee et ah in 
[27] due to the lack of high i polarizations data. Here we derive the distance priors by 
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using the MCMC chains from Planck TT,TE,EE + lowP data in Planck Legacy Archive 
(PLA) [23]. The chains of I a and R can be derived from public MCMC chains by using 
Eqs. Q and ([^. Marginalizing over the remaining parameters, we get the mean values 
and errors of {i?, Ia, kl^lP, Ug] as well as their covariance matrix. 

In Fig. there is a comparison among results derived with the wCDM model but from 
different datasets. We hnd that these three results are consistent with each other, and 
our results (blue ones) impose at least 30% tighter constraints on the distance priors than 
those given by Planck TT+lowP data. 

Meanwhile, in order to explain that the distance priors are effective observables 123. 
we use two different models, i.e. ACDM and tcCDM to derive them respectively. The 
Gaussian likelihood [21] in {i?,/ a, Us} and the covariance matrix are given by the 
upper two sets in Tab. We see that the two sets are stable. 

However, due to the smoothing effect of CMB lensing on the power spectrum [28], 
the third set derived in the curved model (with free Qk) is different from the upper 
two sets by la. Thus, we provide another set of distance priors which are derived with 
free Al (marginalized over later). Here Al denotes the amplitude of the lensing power 
spectrum. This set of distance priors is more conservative and used throughout our 
following discussions. 

^From the PLA, we can find four sets of chains named as 
base_plikHM_TTTEEEdowTEB, base_w_plikHM_TTTEEEdowTEB, 
base_omegak_plikHM_TTTEEE_lowTEB and base_Alens_phkHM_TTTEEE_lowTEB, 
respectively. These four sets of chains are used in our paper to generate Table 1. 
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Figure 1: Marginalized two-dimension probability contours and one-dimension probability 
distribution functions of distance priors derived with the tcCDM model from Planck 2013 
TT-I-WP (green), Planck TT -|- lowP (red) and Planck TT,TE,EE -|- lowP (blue). 


To check our distance priors, we constrain the cosmological parameters in the ACDM 
model by using the distance priors, and compare them with those constrained by the 
global fitting from the Planck TT,TE,EE-|-lowP data. The Xdistance pUors given by 

Xdistance priors ~ ^ ~ di){C )ij{Xj — dj), (8) 

where Xi = {R{z^),Ia{z*), are values predicted by ACDM, di = j^pianck^ 

and Cij is the inverse) is given by the fourth set in Tab. Here we use an ap¬ 

proximate expression of z* to calculate Xi, m, 


z, = 1048|1 + 0.00124(iijV)-"’*®||l + , 

where 

0.0738(Dbh2)-o-238 

“ 1 + 39.5(Dfeh2)0.763 ’ 

0.560 

“ l + 21.1(Dfch2)l-81 ■ 


(9) 


( 10 ) 

( 11 ) 
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ACDM 

Planck TT,TE,EE + lowP 

R 

Ia 

Qbh^ 

ris 

R 

1.7496 ±0.0050 

1.0 

0.49 

-0.69 

-0.77 

Ia 

301.505 ±0.092 

0.49 

1.0 

-0.37 

-0.38 


0.02225 ±0.00016 

-0.69 

-0.37 

1.0 

0.51 

Us 

0.9646 ± 0.0049 

-0.77 

-0.38 

0.51 

1.0 

wCBM 

Planck TT,TE,EE ± lowP 

R 

Ia 

Qbh^ 

Us 

R 

1.7488 ±0.0049 

1.0 

0.49 

-0.68 

-0.78 

Ia 

301.498 ±0.091 

0.49 

1.0 

-0.38 

-0.37 


0.02228 ±0.00016 

-0.68 

-0.38 

1.0 

0.52 

ris 

0.9648 ± 0.0048 

-0.78 

-0.37 

0.52 

1.0 

ACDM+Dfc 

Planck TT,TE,EE ± lowP 

R 

Ia 

Qbh^ 

ris 

R 

1.7449 ±0.0052 

1.0 

0.47 

-0.71 

-0.79 

Ia 

301.465 ±0.093 

0.47 

1.0 

-0.37 

-0.36 


0.02241 ±0.00017 

-0.71 

-0.37 

1.0 

0.54 

Us 

0.9679 ± 0.0048 

-0.79 

-0.36 

0.54 

1.0 

ACDM+Al 

Planck TT,TE,EE ± lowP 

R 

Ia 

Qbh^ 

Us 

R 

1.7448 ±0.0054 

1.0 

0.53 

-0.73 

-0.80 

Ia 

301.460 ± 0.094 

0.53 

1.0 

-0.42 

-0.43 

Qbh^ 

0.02240 ±0.00017 

-0.73 

-0.42 

1.0 

0.59 

Us 

0.9680 ±0.0051 

-0.80 

-0.43 

0.59 

1.0 


Table 1: Distance priors in different cosmological models from Planck TT,TE,EE + lowP. 
Here we also list the scalar spectral index it will be useful if one wants to deal with 
the matter power spectrum. 
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The two constraints on the parameters {^2^, Hq, Qhh^} from the distance priors and the 
global htting are showed in Fig. We hnd that there are almost complete overlaps 
between these two sets. Onr resnlts indicate that the distance priors can be taken as a 
good snbstitnte for the global htting (marginalized over Al) for ACDM model. 



Figure 2: Marginalized constraints on parameters in the ACDM model. Contours given 
by distance priors only are the blue ones, and contours given by Planck TT,TE,EE+lowP 
are the red ones. 


3 Analysis 

3.1 Dark Energy Model 

Although distance priors can be used to constrain the cosmological parameters in ACDM 
model, they can’t constrain the evolving DE model well, such as tcCDM and CPL, because 
DE only plays a crucial role at low redshifts. In order to get a better constraint on the 
equation of state (EOS) of DE, we need to combine some other low redshift observations, 
such as the combination of supernova Union 2.1 compilation of 580 SNe (Union 2.1) [30] 
and low redshift Baryon Acoustic Oscillation (BAO) which are listed in Tab. Here the 













^eff 

measurement 

name 

reference 

0.106 

rs{zd)/Dv = 0.336 ±0.015 

6DFGS 

ra 

0.15 

Dv/rs{zd) = (664 ± 25)/152.66 

MGS 

ra 

0.32 

Dv/r^Zd) = (1264±25)/153.19 

BOSS LOWZ 

ra 

0.57 

Dv/r^Zd) = (2056 ± 20)/153.19 

BOSS GMASS 

ra 


Table 2: The low redshift BAO data. Here the Zd is given by Eisenstein & Hu formula. 


volume-averaged effective distance Dy is dehned by 


Dv{z) = c 


(1 + zYD\{z) 


H{z) 


The baryon drag epoch za is given by Eisenstein & Hu [Mj, i.e. 


Zd = 


1 + 0.659(f2™h2)0'828 




where 


62 = 0.238{nmhy-^^^ . 


( 12 ) 

(13) 

(14) 

(15) 


In this section, we consider a freely binned dark energy (FBDE) model proposed in 
[35113B] . In principle, we can consider a freely N-binned dark energy model. Compared 
to ACDM model, there are more additional free parameters, i.e. the binning redshift 
points and EOS of each bin: Zi, ...,Zi, ...,zn_i and Wi, Unfortunately, it will 

be very time-consuming if all of these parameters are allowed to vary freely in our full 
MCMC analysis. Therefore, we will hrst hnd a set of {zi} to minimize our analysis, 
and then we £x {zi} to do the MCMC sampling of the reduced parameter space. In 
this paper, we deal with a freely 3-binned dark energy model for simplicity. For the hrst 
redshift span {z < zi), the EOS of DE is denoted by wy for the second redshift span 
{zi < z < Z 2 ), the EOS of DE is W 2 . Since at higher redshifts (z > Z 2 ) the energy density 
of DE become quite small compared to that of matter, the expansion of the universe is 
insensitive to DE and the EOS of DE can be hxed as —1 for convenience. Meanwhile, 
we take Z 2 = 1.5, but Zi(> 0.1111) is determined by a best-fit analysis which can yield a 
minimal x^- lu this FBDE model, we have 


Pde(z) 

Pde(O) 


'(l + z)3(i+"'i), z<zi; 

< (l + ^i)3('+"'i) , Z1<Z<Z2; 

/ \ 3(1+1/12) 

^(l + ;^l)3d+-d , ;^>^2. 


(16) 
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Now we can adopt the combined datasets of distance priors, Union 2.1 and low redshift 
BAO to constrain FBDE model with free parameters {VLrn,HQ,wi,W 2 ^^bh‘^-iZi}. We hrst 
set the option “action=2” in CosmoMC to minimize cind we get zi = 0.2256. Then 
we £x Zi = 0.2256 and set the option “action=0” to do a MCMC sampling and we 
constrain the other cosmological parameters. The EOS of FBDE is showed in Fig. 
where wi = —1.0473±0.1852 and W 2 = —0.9720±0.0932 at 68% C.L. which are consistent 
with w = —1 within la C.L.. Moreover, the constraints on and Hq are showed in 
Fig. 1^ Maginallizing over other parameters, we get Ho = 68.34 ± 1.83 km/s ■ Mpc“^ and 
= 0.3043 ± 0.0165 at la C.L.. 



Figure 3: EOS of FBDE in spatially flat universe. 
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Figure 4: Marginalized 2D contours of 1 a and 2 a C.L. in Hq — plane for FBDE given 
by a combination of distance priors, Union 2.1 and low redshift BAO data. 

Here we are also interested in how the former results change in the FBDE model once 
the spatial curvature is taken into account, namely the EBDE+D*, model. Similar to 
the case with Vtk = 0, we get zi = 0.2257 and then £x it to figure out the constraints 
on the other parameters {Q^n, -^o, ^k,'Wi,W 2 , Our results show up in Figs, [^and 

Here wi = —1.1560 ± 0.2178 and W 2 = —0.7977 ± 0.2234 at la C.L., which are also 
consistent with the prediction of ACDM model. The constraints on the other parameters 
are = 0.3017 ± 0.0160, Bq = 68.69 ± 1.80 km/s ■ Mpc-^ and Qk = 0.0053 ± 0.0062 at 
la C.L.. A spatially flat universe is included within 68% C.L.. 
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Figure 5: EOS of FBDE with spatial curvature. 



n,„ H„ fl. 


Figure 6: Marginalized two-dimension probability contours and one-dimension probability 
distribution functions of Hq and for FBDE-|-Ofc model given by a combination of 
distance priors, Union 2.1 and BAO measurements at low redshift. 
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3.2 BAO at z = 2.34 


The BAO feature at high redshift can be measured through absorption in Lya forest 
(LyaF) |37]. According to [38l[39], ^/H{z = 2.34) predicted by the ACDM model con¬ 
strained by Planck 2013 data is approximately 2a lower than that measured by the LyaF 
auto-correlation, while Da{z = 2.36) is approximately 2cr larger than that from the LyaF- 
QSO cross correlation. Furthermore, c/{H{z = 2.34)rs{zd)) = 9.18 ± 0.28 and O.lS^g;^? 
from LyaF auto-correlation and Combined LyaF, respectively. However, the ACDM con¬ 
strained by Planck 2015 data [2H] predicts c/{H{z = 2.34:)rs{zd)) = 8.586 ± 0.021 which 
is more than 2cr C.L. lower than those from both LyaF auto-correlation and Combined 
LyaF. 

Here we wonder whether the tension between Planck 2015 data and the BAO at 
z = 2.34 from LyaF is relaxed in the FBDE and FBDE-|-f2fc models, respectively. The 
two dimensional contours on Da{z = 2.34:)/rs{zd) and c/{H{z = 2.34:)rs{zd)) predicted by 
the FBDE and FBDE-I-Dfc models are given in Fig. where the purple and green crosses 
correspond to LyaF auto-correlation and Combined LyaF, respectively. The marginalized 
Da{z = 2.3A)/rs{zd) and c/{H{z = 2.34:)rs{zd)) are illustrated in Fig. The measure¬ 
ments of LyaF auto-correlation and Combined LyaF and the predictions of FBDE and 
FBDE-|-Dfc models are summarized in Tab. 

To summarize, FBDE and FBDE-|-Dfc models give consistent predictions on both 
Da{z = 2.34l)/rs{zd) and c/{H{z = 2.34:)rs{zd)). But c/{H{z = 2.34l)rs{zd)) predicted 
by FBDE and FBDE-|-Dfc models are still lower than that from LyaF auto-correlation 
measurement at around 2cr C.L., and become even worse compared to that from Com¬ 
bined LyaF. Even though Da{z = 2.34)/rs{zd) predicted by FBDE and FBDE-|-Dfc models 
are consistent with that from Lya auto-correlation measurement, there is still an around 
2a discrepancy compared to Combined LyaF measurement. In a word, relaxing the dark 
energy model from the cosmological constant and adding the spatial curvature cannot 
signihcantly relax the tension between Planck 2015 data and the BAO a.t z = 2.34 from 
LyaF. 
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Figure?: The 68% and 95% constraints on ( 2 : = 2.?>4) /rs{zd) andc/(if(^ = 2.34:)rs{zd)) 
from FBDE (red) and FBDE+fifc (blue). The purple and green crossed denote the mea¬ 
surements from Lya auto-correlation and Combined Lya at la C.L. respectively. 
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Figure 8: Comparison between measurements and our model predictions of c/{H{z = 
2M)rs{zd)) and Da{z = 2M)/rs{zd). 
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Measurement 

Parameter 

mean 

la 

2a 

Lya auto-correlation 

Da{z = 2M)/rs{zd) 

11.28 

+ 0.65 

- 0.65 

+ 2.80 

- 1.20 


c/{H{z = 2M)rs{zd)) 

9.18 

+ 0.28 

- 0.28 

+ 0.60 

- 0.60 

Combined Lya 

Da{z = 2M)/rs{zd) 

10.93 

+ 0.35 

- 0.34 

+ 0.75 

- 0.65 


c/{H{z = 2.M)rs{zd)) 

9.15 

+ 0.20 

- 0.21 

+ 0.40 

- 0.42 

Model 

Parameter 

mean 

la 

2a 

FBDE 

Da{z = 2M)/rs{zd) 

11.72 

+ 0.04 

- 0.04 

+ 0.08 

- 0.08 


c/{H{z = 2M)rs{zd)) 

8.59 

+ 0.03 

- 0.02 

+ 0.05 

- 0.05 

FBDE-|-Dfc 

Da{z = 2M)/rs{zd) 

11.57 

+ 0.23 

- 0.19 

+ 0.40 

- 0.41 


c/{H{z = 2M)rs{zd)) 

8.47 

+ 0.20 

- 0.13 

+ 0.30 

- 0.34 


Table 3: The measurements and predictions of c/{H{z = 2.34)rs(;2d)) and Da{z = 
2M)/r,{zd). 

4 Summary and discussion 

We update the distance priors by utilizing the CMB data, especially including high i 
polarizations data, recently released by the Planck Collaboration. Compared to those 
given in |22] where only Planck TT+LowP data were used, our results impose at least 
30% tighter constraints. Combining the distance priors given in this paper with Union 2.1 
and low redshift BAO datasets, we constrain the cosmological parameters in both FBDE 
and FBDE+f2fc models and hnd that the ACDM model in a spatially flat universe is 
included within Icr C.L.. 

In the literatures, there are debates on the tensions on Da{z = 2.3A)/rs{zd) and 
c/{H{z = 2.?>4:)rs{zd)) between Planck 2015 data and the BAO at = 2.34 from LyaF 
in the ACDM model. We hnd that these tensions cannot be relaxed in both FBDE and 
FBDE+Dfc models. Of course, it is also possible that there are some unknown systemic 
errors in the data analysis for CMB and/or LyaF. It is always worthy investigating all of 
these possibilities in the near future. 
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